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Abstract. Analysis of the multi-wavelength broadband photometry between rest-frame ultraviolet and near- 
infrared wavelengths indicates that the extinction corrected star-formation rates per unit stellar mass of a small 
sample of gamma-ray burst (GRB) host galaxies are higher than those of prototypical, nearby starbursts. This 
result, the confirmed detection of the host of GRB980703 at radio wavelengths, and the tentative evidence in favor 
of a supernova light curve underlying the visible light transient associated with some GRBs provides evidence 
for a connection between stellar phenomenon and GRBs. Fitting population synthesis models to the multiband 
photometry of the host galaxies reveals the presence of a young stellar population with age less than 50 Myr in 4 
out of 6 galaxies, albeit with large uncertainties. Determining the age of the stellar population in a large number 
of GRB host galaxies using this technique could be one of the more reliable ways of distinguishing between the 
collapsar model for GRBs and models that involve the merger of degenerate objects in a binary system. 



THE STARBURST-GRB CONNECTION 

Detailed analysis of the cosmic infrared background and 
long wavelength galaxy counts by various groups have 
illustrated the large (>70%) contribution from infrared 
luminous galaxies (Lir=L(8 — 1000 /jm)>10 n L©) to the 
star-formation rate density at z < 3 [3, 6]. Infrared lu- 
minous galaxies emit as much as 90% of their bolo- 
metric luminosity at far-infrared wavelengths indicating 
that dust reprocessing is extremely significant at z^l — 3. 
Most of these galaxies are inconspicuous at optical/UV 
wavelengths except that many of them show evidence 
of merger activity. This seems to suggest that tidally in- 
duced starbursts dominate the co-moving star-formation 
rate (SFR). However, the classification of a starburst 
galaxy based on its UV-determined SFR is ambiguous. 
For example M82, which is not an infrared luminous 
galaxy, is classified as a starburst galaxy based on stel- 
lar population synthesis model fits to it's optical/near- 
infrared multiband photometry. These fits yield a young 
age for the stellar population, of order 10 Myr. While it's 
star-formation rate is about 7 M Q yr , no higher than 
that of M51 (the Whirlpool Galaxy), it's specific star- 
formation rate i.e. the star-formation rate per unit stellar 
mass is about an order of magnitude higher than M51, 
promoting it to the 'starburst' galaxy category. Thus, if 
GRBs were indeed associated with massive stars, almost 
all the hosts should be starburst galaxies with a large 
fraction of these being infrared luminous objects. To test 
this hypothesis it is necessary to obtain either a measure 
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FIGURE 1. Minimum infrared luminosity that can be de- 
tected by various long wavelength observations as a function 
of redshift assuming the template spectral energy distributions 
of [3]. The sensitivity values adopted correspond to the limits of 
ISOCAM (15,um), SCUBA (850,1ml) and VLA/WSRT (21 cm). 
The 24/jm sensitivity is the MIPS/SIRTF limit for a ~30 minute 
integration. The mid-infrared regime is the most sensitive to 
estimating dust-enshrouded star-formation at z < 2. If we ig- 
nore contributions to the IR light from an AGN, Lir = 10 10 Lq 
corresponds to a star formation rate of 1.7 M0/yr. 



of the far-infrared luminosity or the extinction corrected 
specific SFR in the GRB host galaxies. Figure 1 shows 
the sensitivity of the deepest surveys at a variety of mid- 
and far-infrared wavelengths which allow a direct deter- 



mination of Lir (see also [5]). Clearly, the mid-infrared 
regime (7—25 /jm) is most sensitive to detecting infrared 
luminous galaxies at z < 2.5. However, even the deepest 
surveys at these wavelengths observe only galaxies with 
Lir > 10 11 Lq. 

This is further illustrated by the fact that of the 16 
GRB host galaxies that have been observed so far at 
850 /jm using the SCUBA instrument, only one of them 
(GRBO 10222) has been detected at late-times after the 
afterglow has faded completely [1]. The derived SFR for 
this object is comparable to the SFR inferred from the 
radio observations of GRB980703 [1] suggesting that it 
is an ultraluminous infrared galaxy (ULIG). 

The low detection rate of GRB host galaxies in the 
submillimeter questions the validity of the association 
between star-formation and GRBs. However, at the typ- 
ical SCUBA sensitivity of 3o ^3 mJy, only the bright- 
est ULIGs that are more luminous than L iR ^4x 10 12 L Q 
can be detected at z > 1 . These hyper luminous galaxies 
comprise only a small subset of the entire ULIG pop- 
ulation. Since the total contribution from ULIGs to the 
global star-formation history is at most 30% [3] the con- 
tribution from these hyper-luminous sources is less than 
10%. 

Thus, as inferred independently by Ramirez-Ruiz et 
al. [10], only about 10% of GRB hosts are likely to be 
detected by SCUBA at sensitivity limits of 3 mJy which 
is consistent with the small detection rate mentioned 
above and the small number of statistics. MIPS/SIRTF 
observations will be vital in deriving the far-infrared 
luminosity of a large number of GRB hosts. 

Fortunately, at the faint end of the infrared luminos- 
ity function, a complementary technique which utilizes 
the UV-slope seems to measure the dust-obscured star- 
formation reasonably well [9]. At Lir < 4 x 10 11 L , 
much of the dust obscuration is optically thin i.e. the far- 
infrared emission is dust reprocessed UV emission from 
the same part of the galaxy. However, in the more lu- 
minous sources the dust is concentrated so strongly that 
the UV opacity is very high in the neighborhood of the 
star-forming regions and so all the UV photons are re- 
radiated at longer wavelengths. In this case, the UV and 
far-infrared observations trace the low and high opacity 
regions of the same galaxy respectively and the UV-slope 
is only a lower limit to the amount of obscuration. 

To search for signatures of a starburst in the GRB 
hosts using this technique, the multiband photometry on 
GRB host galaxies between rest-frame ultraviolet and 
near-infrared wavelengths were fit using the stellar pop- 
ulation synthesis models of Bruzual & Chariot, taking 
into account internal dust extinction [4]. The total (ob- 
scured+UV) star-formation rate in these galaxies is not 
found to be unusually high (Figure 2). Although only 2 
of the 6 host galaxies with accurate photometry have de- 
rived Lir > 10 11 L© and high resultant SFRs, the ratio 
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FIGURE 2. Total (obscured+unobscured) star-formation 
rates for selected gamma-ray burst host galaxies derived using 
the )3-slope technique plotted against the ratio between the star 
formation rates and the midpoint of the range of stellar masses 
derived from fits to the multiband photometry. The uncertainty 
on the GRB971214 point is quite large because of the lack of 
good photometry at multiple wavelengths. The star-formation 
rates for the GRB hosts are lower limits, so the data points are 
likely to move higher and to the right. Estimates of the stellar 
mass typically have 95% confidence intervals that span an order 
of magnitude. Also plotted are the corresponding values for two 
prototypical starbursts Arp220 and M82 and the relatively qui- 
escent Sbc galaxy M5 1 . GRB hosts have star-formation rates 
per unit stellar mass much higher than local starbursts. 



between the dust-obscured SFR and the unobscured SFR 
for the GRB hosts has an average value of ^4. This is in 
agreement with the redshift dependent value of 3— 7 that 
Chary & Elbaz [3] find for the ratio between the global 
comoving dust-obscured star-formation rate and the un- 
obscured star-formation rate derived from the UV. 

The near-infrared wavelengths constrain the stellar 
mass and age of the galaxy very well (Figure 3). Our 
analysis of the stellar population age and baryonic mass 
of 6 of the host galaxies with accurate photometry 
and redshifts reveal that all of them have high specific 
star-formation rates compared to local starbursts (Fig- 
ure 2). The spectral energy distribution of 4 of them 
(GRB970228, GRB970508, GRB980613, GRB980703) 
are dominated by a young (age<50 Myr) stellar popu- 
lation which favors the collapsar model. However, 2 of 
them (GRB990123, GRB991208) have stellar popula- 
tions of age ^200 Myr suggesting that the merger of dou- 
ble degenerate systems might power some of the GRBs. 



FUTURE WORK 

Various techniques have been suggested to constrain 
the progenitors of GRBs: studying the environments of 
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FIGURE 3. Optical/near-infrared photometry for the host 
of GRB980703 along with a template SED of a galaxy with 
an exponentially decaying starburst. The template consists of 
a young stellar population fit to the optical/UV photometry 
and the maximal contribution from an older stellar population 
which might have formed in an initial burst at z ~10 and 
has evolved over the 5 Gyr between z = 10 and z = 1.097. 
The solid black line is the sum of the two components. The 
luminosity of many GRB host galaxies are dominated by a 
young (10-50 Myr) stellar component which favors an origin 
of GRBs from collapsars. GRB host galaxies whose SEDs are 
dominated by an older stellar component would favor alternate 
mechanisms such as merging neutron stars. The lower plot 
shows the range of values for the extinction and age which 
result in SEDs that are within 2a of the data points at B,V,R,I 
and K. The different contours illustrate the variation in the 
parameters with photometric uncertainty. 



GRBs through spectral analysis of the X-ray afterglow 
[8], searching for signatures in the light curve of the opti- 
cal transient associated with the GRB [2, 7] and studying 
the nature of the host galaxies. Of these, the latter is less 
direct but relatively simple and robust. 

Associating GRBs with stellar phenomenon necessi- 
tates that the host galaxies are either undergoing active 
star-formation or have been through epochs of substan- 
tial star-formation in the past. Since much of the high- 
redshift star-formation is obscured by dust, accurately 
measuring the SFR requires tracing the amount of inter- 
nal extinction in the host galaxies. 

Previous star-formation episodes can be constrained 
either by measuring the metallicity of the host galaxies 
or by determining the stellar mass. Estimates of galaxy 
masses derived by fitting template SEDs to the photom- 
etry between UV and visible wavelengths only provide 
lower limits since they trace the young stellar component 



which has a low mass to light ratio. Accurate masses can 
only be derived by including data at rest-frame X ~ 2 /urn 
which traces the old/cool stellar component. At high 
redshift, IRAC/SIRTF observations will be able to con- 
strain the contribution from this component. Thus, by fit- 
ting population synthesis models to the rest-frame UV to 
near-infrared light one can estimate the mass fraction of 
gas that has gone through the most recent starburst, the 
internal extinction, the age of the starburst and fraction 
of galaxy mass incorporated in an old stellar population. 

Determining the distribution of stellar ages in a 
large number of host galaxies will provide one of 
the more reliable ways to distinguish between models 
where massive stars evolve to become GRBs. In par- 
ticular models which involve the production of GRBs 
through the core-collapse of isolated massive stars into 
a black hole/accretion disk system would take place on 
timescales comparable to the epoch of star-formation 
i.e.^lO Myr. On the other hand merging double degen- 
erate models suffer an evolutionary delay from the epoch 
of star-formation which is induced by the timescale for 
orbital decay of the binary system >100 Myr. 

Thus, if a majority of the GRB hosts are found to 
have high specific SFRs (SFR/stellar mass) and young 
(age<10 Myr) stellar populations, this would be con- 
vincing evidence for 'collapsars' to be the progenitors 
of GRBs. However, most of the high redshift infrared lu- 
minous galaxies show the presence of a significant older 
stellar population that is dominating the rest-frame near- 
infrared light of the galaxy. This indicates that infrared 
luminous galaxies have undergone previous epochs of 
star-formation which could potentially provide an origin 
for double degenerate systems (Figure 3). As a result, if 
the majority of the hosts have a large old stellar popu- 
lation in addition to the young component, this would 
make it impossible to discriminate between the two mod- 
els using this technique. A positional analysis of the burst 
locations on the host would then have to used as the dis- 
criminator due to the large kick velocities imparted to the 
degenerate objects from their supernovae. 
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